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We study the evolution of entanglement of a pair of coupled, non-resonant harmonic oscillators 
in contact with an environment. For both the cases of a common bath and of two separate baths 
for each of the oscillators, a full master equation is provided without rotating wave approximation. 
This allows us to characterize the entanglement dynamics as a function of the diversity between 
the oscillators frequencies and their mutual coupling. Also the correlation between the occupation 
numbers is considered to explore the degree of quantumness of the system. The singular effect of 
the resonance condition (identical oscillators) and its relationship with the possibility of preserv- 
ing asymptotic entanglement are discussed. The importance of the bath's memory properties is 
investigated by comparing Markovian and non-Markovian evolutions. 

PACS numbers: 03.65.Yz, 03.65.Ud 



I. INTRODUCTION 

Coupled harmonic oscillators are the first approxima- 
tion to a broad class of extended systems not only in 
different fields of physics but also in chemistry and biol- 
ogy. Within the quantum formalism, they are the basis 
of the description of electromagnetic field interactions in 
quantum optics and approximate lattice systems in dif- 
ferent traps in atomic physics [3, Hi- Moreover, in the last 
few years there has been impressive progress towards the 
cooling and back-action evasion measurement of 'macro- 
scopic' -in terms of number of atoms- harmonic oscilla- 
tors, allowing the observation of their quantum behavior. 
Two main class of systems experimentally realized are na- 
noelectromechanical structures (NEMS) [!, 0] and differ- 
ent kinds of optomechanical systems (BJ where nano- and 
micromechanical devices, cavities or suspended mirrors 
are respectively coupled to single electrons or light. As an 
example, it has recently been reported the observation of 
NEMS extremely near to the ground state of motion with 
an occupation factor of just 3.8 [|[. These experiments 
would allow to observe coherent quantum superposition 
states, entanglement and to study decoherence processes 
in a controllable way on massive solid-state objects. 

Phenomena associated to the coupling of these quan- 
tum oscillators have been revisited within the context of 
quantum information in many theoretical studies during 
the last decade [rj-[T2l|. Entropy and entanglement in ex- 
tended systems with many degrees of freedom (harmonic 
chains or lattices) @, have been characterized in funda- 
mental and thermal states, exploring scaling laws and 
connections with phase transitions Q ■ An important ad- 
vantage is that these systems admit Gaussian state solu- 
tions having a well defined 13[ computable [14| measure 
of entanglement, the logarithmic negativity [15|. The 
question of the generation of entanglement has also been 
addressed considering oscillators whose parameters are 
modulated in time [§, 9] . In these studies losses were gen- 
erally neglected, while recently the effects of decoherence 
on a pair of entangled oscillators have been considered in 



presence of dissipation through baths of infinite oscilla- 
tors (Tol - fl2| . Our aim in this paper is to analyze a rather 
unexplored aspect of this problem that is the effect of 
the diversity on the entanglement between coupled har- 
monic oscillators in different situations. Indeed, instead 
of considering identical oscillators, we look at the effects 
of detuning between their frequencies, lj\ and ui2- 

The interest about diversity effects on entanglement is 
both theoretical and related to experimental issues. We 
mention, for instance, the effect of diversity of frequen- 
cies of two photons entering in a beam splitter. This 
leads to a completely different output with respect to the 
case of indistinguishable photons [III [I?} [H( Coupling 
between different harmonic modes has been also exten- 
sively studied in quantum optics in presence of nonlin- 
ear interactions and parametric coupling, allowing, for 
instance, generation of entan glem ent between photons 
pairs and intense light beams (loL [20| of different colors. 
In that context, however, the frequency diversity of each 
pair of oscillators is compensated by a third mode and 
is in general not relevant, while here we focus on pairs 
of mechanical oscillators off-resonance and with constant 
couplings. The main expected consequence is indeed an 
effective decoupling of the oscillators due to fast rotation 
of their interaction term and we will show the effects on 
the robustness of entanglement. 

It is clear that the identity of the oscillators in general 
is a very peculiar and strong assumption not always jus- 
tified and introducing a symmetry into the system with 
deep consequences. It is indeed important to clarify the 
effects of relaxation of this symmetry introducing some 
diversity. In an extension from two coupled oscillators to 
an array, this will imply a break in the translational sym- 
metry. Apart from the fundamental interest, we point 
out that in experimental realization of engineered arrays 
of massive quantum oscillators, some diversity between 
them might actually be unavoidable. Coupled oscillators 
with different frequencies have been also suggested for 
quantum limited measurements (2lj . Finally, in many 
experiments coupled oscillators actually model different 
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physical entities (for instance radiation and a moving 
mirror in optomechanics) and a symmetric Hamiltonian 
would describe only a very special case [22^ . 

Our analysis of entanglement evolution encompasses 
both diversity between harmonic oscillators and dissipa- 
tion. Once the oscillators are prepared in some entan- 
gled state (for instance through sudden switch of their 
coupling), we look at its robustness increasing diversity 
and coupling strength. It is well known that an object 
whose all degrees of freedom are coupled to an environ- 
ment will decohere into a thermal state with the same 
temperature as the heat bath 0]. The thermal state, un- 
less temperatures arc very low, is separable and highly 
cntropic, so that after thermalization all entanglement 
shared between the coupled oscillators disappears (23j | . 
As the dissipating oscillators reach a separable state in 
a finite time, and not asymptotically, the name of sud- 
den death has been suggested [24 |. Still under certain 
conditions on the ways in which the oscillators dissipate, 
entanglement can survive asymptotically pBj . The tran- 
sition from a quantum to a classical behavior in these 
systems does not convey only fundamental interest being 
also important in view of applications of harmonic sys- 
tems operating at the quantum limit, including quantum 
information processing and not shot-noise limited mea- 
surements of displacement, forces or charges [26j]. These 
phenomena have generally been studied for identical os- 
cillators and, in particular, previous works provide a mas- 
ter equation description in the case of a fully symmetric 
Hamiltonian [l(| E[ ■ Off-resonance oscillators in pres- 
ence of a common bath have recently been considered in 
Ref. [l2j showing that in this case lower temperatures 
are needed to maintain entanglement asymptotically. In 
this work we analyze systematically the role of diver- 
sity on the entanglement dynamics and provide the full 
master equations for w\ ^ UJ2 (i) both for common and 
separate baths for the two oscillators, (ii) without the ap- 
proximation of rotating waves, generally assumed when 
modeling light fields, and (iii) comparing results to the 
non-Markovian case. 

Recently there have been some works focusing on mem- 
ory effects and non-Markovianity. In the case of contin- 
uous variables, non-Markovian effects on entanglement 
evolution have been discussed [HI H3|- ^ n ^ e ^- P^t - 
the authors considered two identical oscillators coupled 
to separate baths in the very high temperature regime 
and analyze how the matching between the frequency of 
the oscillators and the spectral density of the bath af- 
fects both Markovian and non-Markovian dynamics. In 
Ref. [l0| . the non-Markovian evolution has been com- 
pared with the case where both the Markovian limit and 
the rotating wave approximation in the system-bath cou- 
pling are taken. Since the effects of these two approxi- 
mations are not easily separable, it is difficult to single 
out non-Markovian corrections from that analysis. Our 
study of entanglement dynamics in presence of diversity 
is also extended, for the sake of comparison, to the non- 
Markovian case, looking at both common and separate 



baths and showing that deviations are actually negligible. 

In Sect. |n] we introduce the model of a pair of os- 
cillators with different frequencies and coupled through 
position both between them and with the baths of oscil- 
lators. Both master equations for one common and two 
separate baths arc presented, providing all the details in 
the Appendix. Temporal decay of entanglement between 
the oscillators and correlations between the occupation 
numbers are shown in Sect. IIIII analysing the role of 
frequency diversity when also their coupling strength is 
varied. Apart from the entanglement dynamics, we also 
discuss its robustness (asymptotic entanglement for com- 
mon bath) in the context of the symmetry of the system, 
Sect. IIII Bl Non-Markovian deviations from these re- 
sults are shown in Sect. IIVI and further discussion and 
conclusions are left for Sect. [V] 



II. MODEL AND MASTER EQUATIONS 

We consider two harmonic oscillators with the same 
mass and different frequencies coupled to a thermal bath. 
As discussed in Ref.[2g, depending on the distance be- 
tween the two oscillators, different modclizations of the 
system-bath interaction can be done. We will discuss the 
case of two distant objects, which amounts to considering 
the coupling with two independent baths, and compare 
it with the zero-distance scenario (common bath). An- 
alyzing the role of diversity on the quantum features of 
this system, we also generalize previous works on identi- 
cal oscillators dissipating in common and separate baths 

A. Separate baths 

The model Hamiltonian, with each oscillator coupled 
to an infinite number of oscillators (separate baths), is 
H se P = h s + H s B ep + H s s e g. The system Hamiltonian 

rr Pi , 1 2„2 , Vl , 1 2„2 , \ /1\ 

3 = y 2 UJiXi y 2 W2 ' T2 * 2 ^ ' 

describes two oscillators with different frequencies oji.2, 
coupled through their positions, 

2 /p(i)2 , \ 

ff B p = EE V + (2) 

is the free Hamiltonian of two (identical) bosonic baths, 
and 

^f = E4 1) 4V+E4 2) 4 2) -2 (3) 

encompasses the system-bath interaction. 

The master equation for the reduced density matrix of 
the two oscillators, up to the second order in H S gg (weak 
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coupling limit) is, assuming h = 1, 



dp 
If 



Ttr 1 



1 2 



(4) 



In Appendix we give an explicit derivation of the master 
equation together with the definition of the coefficients. 
Then, p is subject to energy renormalization (e,j), dis- 
sipation Tij , and diffusion Dij , Fij . Coefficients in Eq. 
(j4]) depend on time, but in the following, we will con- 
sider the Markovian limit (performed sending t — > oo 
in Eqs. (|A17IIA20)) ). Non-Markovian corrections will be 
discussed in Sec. IIVI We anticipate that we will focus on 
the Markovian limit because corrections dropping from 
such an approximation turn out to be negligible in most 
of the cases, both for equal or different frequencies oji 
and u>2; and for common or separate baths. 

To obtain an explicit expression for ey , Dij , Tij , Fij , we 
need to know the density of states of the baths, defined 
for both of them as 



J(Q) 



\ 2 

k K 



<s(n-n fc ). 



(5) 



Here, we will consider explicitly the Ohmic environment 
with a Lorentz-Drude cut-off function, whose spectral 
density is 



A 2 



ft 2 



(0) 



We learned from the master equation (j4j that bare fre- 
quencies of the oscillators are renormalized because of 
the presence of eij, and this renormalization turns out 
to depend on the frequency cut-off A. This undesirable 
unphysical effect can be removed by adding to the initial 
Hamiltonian counter-terms [29j which exactly compen- 
sate the asymptotic values of ey. This is accomplished 
by replacing Cij(t) with Cij(t) — (oo) in Eq. Q. Then, 
in the Markovian case, we will simply drop them from 
the master equation. This renormalization procedure 
amounts to redefining the natural frequency in terms of 
the observed frequency. 



B. Common bath 



We introduce now the case of common bath modeled 



by 



2 2 



(7) 



and we consider the same spectral density ((5]). The in- 
teraction term reads 



H SB = X] ^Xk(xi + x 2 ), 



(8) 



meaning that the system is coupled to the bath through 
the mode x+ = (x\ + X2)/v2, while the mode X- 

i — X2)/V% is decoupled. The master equation in this 
case is 



dp 
~dt 



-i[H s ,p\ - - ^(4 - <].■■'■ 

1 2 

--j=^2{ie%[xiX + ,p] + D u [x + , [xi,p]] 
+iT u [x + ,{pi,p}] - F u [x+, \pi,p]]}, (9) 



(see the Appendix for the definition of coefficients) . Cou- 
pled systems dissipating in a common bath have been 
recently subject of interest because of the possibility to 
preserve entanglement asymptotically at high tempera- 
tures, being this a major distinctive feature with respect 
to the case of separate baths. In Sect. IIII B[ we will 
study entanglement evolution under Eq. ((9]) revisiting 
the possibility to maintain entanglement asymptotically 
in presence of diversity (frequency detuning) . As a mat- 
ter of fact, whereas in the case of equal frequencies the 
decoupled mode X- is an eigenmode of the isolated sys- 
tem, this is no longer true once lo\ ^ ui2- We will see how 
asymptotic entanglement can be preserved in presence 
of diversity, through a proper choice of the system-bath 
coupling constants. 

Another difference in the case of common bath is the 
lack of the symmetry A —> —A. A closer look to the form 
of the total Hamiltonian in presence of separate baths 
allows us to single out its symmetry properties. H sep is 
clearly modified by the transformation A — > —A. How- 
ever, once the canonical transformation U is introduced, 
such that Ux2U^ = —x 2 and Up2U^ = —p2, it is easy 
to show that the master equation is invariant under the 
combined action of U and the flip of A. In fact, while Hs 
is left unmodified, the change in Hg e £ (that would cor- 

(2) (2\ 

respond to flip X k in —X k ) is ineffective with respect 
to the evolution of the reduced density matrix, being p 
determined only by contributions proportional to (A^) 2 . 
A simple consequence of this symmetry can be observed 
studying the evolution of the two-mode squeezed state 



\*tms}= v/l - T^p n/2 \n)\n), 



(10) 



n=Q 



where /i = tanh r, and r is the squeezing amplitude. 
For this state, the canonical transformation U amounts 
to changing r in — r. Then, the same evolution must be 
obtained considering A in H sep and a given squeezing r 
in the initial condition, or —A and — r. All these con- 
siderations are true only in the case of separate baths. 
Once the common bath is taken, because of the action 
of U, the mode coupled to the bath is x_ instead of x+. 
Therefore, in the case of separate baths, results for oppo- 
site A are equivalent to a change in the sign of the initial 
squeezing, while this is not the case for a common bath. 



4 



III. ENTANGLEMENT AND QUANTUM 
CORRELATIONS 

The calculation of bipartite entanglement for mixed 
states is generally an unsolved task. Nevertheless, the 
criterion of positivity of the partial transposed density 
matrix of Gaussian two-mode states o Tb is necessary and 
sufficient for their separability [l]J[3(l|3l|. The amount of 
entanglement can be measured through the logarithmic 
negativity, defined as Ejsf = log 2 ||p Tf5 ||, being ||p Tf5 || the 
modulus of the sum of the negative eigenvalues of p Tn 
For any two-mode Gaussian state, the logarithmic 
negativity is E^f = max[0, — log 2 2A_], where A_ is the 
smallest symplectic eigenvalue of p Tn . 

Because the system is connected to heat baths, the 
dissipating degree of freedom will ultimately attain the 
thermal state corresponding to the system Hamiltonian. 
The properties of entanglement in the thermal state are 
well studied, whence it is known that only for very low 
temperatures entanglement is present between the two 
oscillators. Therefore we focus our attention on the time 
evolution of the decohcrence process itself, namely how 
fast an initially entangled state decays into a separable 
one, for separate (Sect. IIII A|) and common (Sect. MI B[) 
baths. 



A. Entanglement decay time 

We consider as a non-separable initial state the two- 
mode squeezed vacuum (|10[) . with squeezing parame- 
ter r. This state has an amount of entanglement pro- 
portional to r and consists in the orthogonal modes 
x± = (xi ± cc2)/v / 2 being simultaneously squeezed and 
stretched, respectively, by an amount r. Such a state, for 
infinite squeezing, gives the maximally entangled state 
for continuous variables, long known due to the famous 
paper by Einstein, Rosen and Podolsky [32j . 

The entanglement evolution in presence of separate 
baths is shown in Fig. Q}i for different parameters choices, 
showing the effects of increasing the coupling strength 
and the detuning. Entanglement, in general, decays with 
an oscillatory behavior with oscillation amplitude de- 
pending on the A strength Q. We study the last time 
at which entanglement vanishes, tp. In Fig. [5] and [3] we 
scan how long it takes for such an entangled state (with 
r = 2) to become separable when the oscillators detun- 
ing and coupling strength are varied. Indeed, the last 
time tp is represented in Fig. [HJfor different values of the 
ratio between frequencies {UI2/001) and coupling between 
oscillators (A/cof). 

For this choice of parameters, Fig. [5] shows that en- 
tanglement is present for few tens of periods. Due to the 
scaling of temperature, time, and couplings with the fre- 
quency of the first oscillator, cji, it is clear that results 
shown in Figs. [5] and [3] arc not symmetric by exchange 
of the oscillators role. In other words, even if it is physi- 
cally equivalent to have the first oscillator with half the 



frequency of the second one or the second half the first 
one, due to the scaling, these figures are not symmetric 
with respect to the line 102/ oo\ = 1. 

Figure [5] shows two main features: (i) an increase of 
the survival time tF when 102/001 is increased; (ii) an in- 
crease of tp with |A| in absence of diversity (co2 = <^i), 
while for 002 > 2wi this dependence is weakened. This 
can also be appreciated in the perspective presented in 
Fig. EI Both features can be explained in terms of the 
eigenmodes of the system Q± (|A10|) and their eigenfrc- 
quencies £l± (|A9[) . Since the eigenmodes do not interact 
with each other, they can be regarded as independent 
channels for decoherence. The eigenmodes will only in- 
teract with near-resonant frequencies in the bath. In ad- 
dition, all variances at the thermal states they approach 
are dependent upon the fraction Q±/2kBT. 

As far as it concerns the feature (i), in absence of cou- 
pling (A = 0, implying = LO2.1) and for L02 — > 00, the 
'effective temperature' of the final thermal state reached 
by the eigenmode Q+ will vanish, T e g^ + = ksT '/0J2 — > 0. 

The eigenmodes will therefore reach respectively a 
thermal state (T c ff.._ = hsT/coi) and a ground state 
(T ff. j+ = 0). Since the presence of entanglement is a 
competition between the reduced purities of individual 
oscillators and the total purity [33| . just by improving 
the final purity of oscillator Q + its time evolution has an 
overall higher purity, thus making entanglement higher 
(and its survival time longer). This effect is equivalent 
to reduce the real temperatures of the baths (while keep- 
ing everything else constant). While the decoherence is 
given by the coupling 7 and is hence not reduced, the 
reached final state is purer and thus entanglement is seen 
to survive longer. 

As far as (ii) is concerned, for L02 = ^1 the eigenfrc- 
quencies of the system are fi± = \/ <-o\ ± |A|. By in- 
creasing |A| up to uo\ the eigenfrequency f2_ vanishes 
and the amount of bath modes which have a similar 
frequency, given by the spectral density, also vanishes 
( J(f2_ —>•())—>• 0). That is, the amount of bath modes 
with which this degree of freedom interacts tends to zero, 
and therefore the amount of decoherence suffered. This 
effect, however, is partially compensated by the fact that 
a vanishing f2_ would imply that this mode will reach a 
thermal state of effective infinite temperature. Were it 
not so, the entanglement would survive asymptotically, 
as in the case of common bath, where for L02 = u>i the 
mode tt- docs not decohere. 

We then see that all major effects related to entan- 
glement decay between coupled oscillators in presence of 
heat baths can be explained in terms of: 1) eigenfrc- 
quencics (where they lie within the spectral density of 
the heat baths), and 2) effective temperatures reached 
by the eigenmodes after thermalization. 

A minor feature is that the increase of tp with A and 
L02 is not completely smooth due to entanglement oscilla- 
tions and this gives rise to the 'arena'-shaped dependence 
that is clearer in Fig. [3j In general, the dependence of 
the decoherence time on the oscillation frequency of the 
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system is negligible for optics experiments at ambient 
temperatures {T e g ± = kgT/£l± ~ 0), but becomes ex- 
tremely important in presence of the lower frequencies of 
mechanical oscillators [34j . 

Notice that we have restricted the coupling within the 
boundary |A| < U1LO2 (triangle in the lower part of Fig. 
[2]), otherwise one of the eigenfrequencies (fi-) would be- 
come imaginary, meaning that trajectories would be un- 
bounded which, though not unphysical in principle, leads 
to leaks in any kind of experiment. 
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FIG. 1: (Color online). Time evolution of an entangled two- 
mode squeezed initial state (|10[) with r = 2, in the case of 
a) separate baths and b) common bath. The parameters are 
ksT — IOcji, 7 = O.OOIcji, and cutoff frequency A = 50o>i. 
The four curves correspond to the points in figures [2] and [4] 
namely A(dotted black), B(solid black), C(dotted orange) and 
D (solid orange). Two features are apparent: in presence of a 
nonzero coupling A (continuous lines) entanglement oscillates 
strongly, and second, in the case of separate baths any initial 
entanglement vanishes fast, while in the case of common bath 
it can be made to survive by having identical frequencies. 
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FIG. 2: (Color online). Parameter scan of £f (see text) as a 
function of the oscillators coupling and frequency detuning, in 
the case of separate baths. We have restricted the scan within 
a limited detuning (ui2 < 2wi) as a result of the absence of 
any important features outside this region. In addition, the 
condition A < ensures reality of the eigenfrequencies in 

the problem. The behavior of If is basically monotonic in 
UJ2/0J1, with a superimposed arena-like shape, coming from 
the oscillatory nature of entanglement. See text for more 
details. The dots A, B, C, D are given in figure[T]as examples. 




B. Common bath 

Let us compare now the previous results with the os- 
cillators evolution with dissipation through a common 
bath, Eq. ((9]). This model can be considered a limit 
case while separate baths would model the opposite one, 
when baths are completely uncorrelated. The presence 
of asymptotic entanglement in this situation has been 
studied by [H O > 

showing that entanglement can 
survive for sufficiently low temperatures (or conversely 



FIG. 3: (Color online). Three dimensional representation of 
figure [2] Here the arena- like shape can be better appreciated. 

for high enough initial squcezings). In the case with dif- 
ferent frequencies Paz and Roncaglia noticed in Ref. [l2| 
that this asymptotic behavior disappears as the detuning 
is increased, meaning that it is highly dependent on the 
frequency matching between oscillators and the fact that 
each oscillator's bath can be regarded as 'perfectly cor- 
related' to the other one. In physical realizations these 
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FIG. 4: (Color online). Parameter scan of tp (see text) as a 
function of the oscillators coupling and frequency detuning, in 
the case of common bath. As seen here, tp is huge whenever 
the detuning is small, and in the order of tens of periods for 
high detuning. As expected, at resonance one of the eigen- 
modes decouples from the baths, leading to an infinite £f- We 
have truncated the plot at heights of oj\tp = 210, otherwise 
this 'mountain riff' is infinitely high (but only strictly infinite 
when cui = LU2). 
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FIG. 5: (Color online). 3D representation of Figure[4] Apart 
from the 'riff' at equal frequencies, the shape is very similar to 
the case of separate baths, with an arena-like bowl shape. We 
have drawn a mesh to guide the reader's eye. Again, heights 
over tp = 210/u\ have been truncated (yellow region). 



two conditions will hardly be met and it is interesting 
to know the effect of deviations from it. We have then 
studied the effects of frequency diversity in the case of 
a common bath, looking again at the robustness of en- 
tanglement in terms of the decay time tp (Figlljj) . In 
Fig. [5] we show results equivalent to what represented in 
Fig. [3] but for common bath. There we see that in the 
resonant case entanglement never vanishes, and so the 
survival time is in that case infinite. In Fig. [5] diverging 
times tp are recognized along the line u\ = uj^. 

When the oscillators frequencies begin to differ, sur- 
vival diminishes fast to values similar to the uncorrclatcd 
baths case. This phenomenon can be understood from 
the master equation j9]). The fact that the coupling to 
the bath is 7(0:1 + X2) establishes that the mode x_ re- 
mains decoupled at all times, thus keeping its coherence 
and contributing positively to a nonzero entanglement. 
If the frequencies arc different, that mode gets coupled 
to x + which is itself coupled to the bath, and its deco- 
herence is transferred to x_. This way, both modes are 
decohered and end up in a thermal state, with no pres- 
ence of entanglement. Hence the survival time is finite. 
Because of this, it is worthy to stress the fact that even if 
the frequency difference is infinitesimal, decoherence will 
eventually show up, even though it does so in a time in- 
versely proportional to the infitcsimum. Thus, increasing 
the difference in frequencies just exacerbates the transfer 
velocity of decoherence from mode x+ to mode X- . This 
argument evidences how artificial the 'equal frequencies' 
assumption is. 

We point out that asymptotic entanglement can be 
found also in presence of frequency diversity and is not 
related to the symmetry present for uj\ = L02 ■ As a matter 
of fact, if the couplings to the common bath are different 



for each oscillator, i.e. 71 7^ 72, we can restore asymp- 
totic entanglement even off-resonance, for lji 7^ 012- By 
noticing the relation between X12 with the eigenmodes 
of the system Q±, we see that the mode Q+ can be un- 
coupled from the bath when the angle 8 fulfills 

cos6» = ^J^= (11) 

with 9 = 9(lui,uj2, A) given in the Appendix. Thus, half 
of the system is kept in a pure state, and we can show that 
entanglement survives asymptotically exactly in the same 
fashion as it did for equal frequencies above. In other 
words, with different couplings to the baths together with 
different frequencies, the system can retain asymptotic 
entanglement if Eq. (fTl) is fulfilled. In any case this phe- 
nomenon is as unique as the 'equal frequencies' one, and 
requires fine-tuned parameters (tuning bath couplings or 
frequency difference) in order to be observable. It should 
then be regarded as exceptional too. Our argument clar- 
ifies that asymptotic entanglement is not a consequence 
of a symmetry in the system, but comes from finely tuned 
decoherent-free degree of freedom of the system. 

C. Twin oscillators 

We have considered entanglement as an indicator of 
the quantumness of our system, measurable for the fam- 
ily of (Gaussian) states considered here. The discrimi- 
nation between the predictions of classical and quantum 
theories for coupled harmonic oscillators has been largely 
studied also in optics. In that context there have been 
many theoretical predictions experimentally confirmed, 
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focusing on the violation of different classical inequali- 
ties, or the positivity of variances [35j |. An example is 
the variance of the difference of the occupation numbers: 



(m 



n 2 ) 



(12) 



where, as usual, rii is the occupation number opera- 
tor of each oscillator, d has been considered in two 
mode squeezed state generated by parametric oscillators 
in optics, to characterize twin beams The quan- 

tum character of the correlations in the occupation num- 
bers comes from the negativity of the variance d and 
can only follow from the negativity of the corresponding 
quasi-probability, the Glauber-Sudarshan representation 
in this case. In analogy with the optical case, we consider 
here twin oscillators looking at the temporal dynamics 
of d. The correspondent fourth order moments can be 
obtained from the covariance matrix, because the states 
we are dealing with are Gaussian. In Fig. [6] we plot 
two examples comparing the evolution of this variance 
and entanglement for the common/separate baths cases. 
We have taken only the negative part of this correlation 
-identifying quantum behavior- and inverted it (plotting 
max(0, —d) ) for ease of comparison with entanglement. 
For the initial entangled state (|10[) . d = —2^/(1 — fj,), 
always negative. In other words, for squeezed states the 
occupation number of one oscillator determines the other 
one. We observe that starting from this value d decays 
with large oscillations and that the sudden deaths of en- 
tanglement and of this correlation coincide up to a frac- 
tion of a period. Still entanglement evolution is smoother 
and twin oscillators temporarily loose their quantum cor- 
relations even for entangled oscillators. 
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FIG. 6: (Color online). Comparison between the time evo- 
lutions of entanglement (black) and the correlation measure 
max(0, — (: (ni — 112) 2 :)) (orange), for the case of a) separate 
baths and b) common bath. We have used equal frequen- 
cies, coupling A = O.lui, damping 7 = O.OOIcji, temperature 
UbT /u)\ — 10 and initial squeezing parameter r = 0.5. The 
envelope of the correlation measure seems to be closely related 
to the entanglement, except for a delay of half a period. 



IV. NON-MARKOVIAN EVOLUTION 

Up to now we have considered the Markovian master 
equation, in the sense that the time-dependent coeffi- 
cients related to the heat bath have been replaced by 
their asymptotic value, obtained by integrating up to in- 
finite time. This implies a complete absence of memory 
in the bath, which does not retain instantaneous infor- 
mation on the dynamics of the two oscillators. However, 
this simplification can be dropped, and we are left with 
a complete description within the weak coupling approxi- 
mation. In this section we analyze how these corrections 
affect the entanglement evolution we discussed so far. 

In Ref. |27j . a comparison has been made consider- 
ing temperatures two orders of magnitude greater than 
the frequency cut-off. The authors found that, when the 
frequency of the oscillators falls inside the bath spectral 
density, entanglement persists for a longer time than in 
a Markovian channel. When there is no resonance be- 
tween reservoir and oscillators, non-Markovian correla- 
tions accelerate decoherence and generate entanglement 
oscillations. 

As it is known, when all the important frequencies 
are much lower than the cut-off, the density of states 



of Eq. © is expected to generate almost memorylcss 
friction |36| . Then, taking u>\ <C A, we compared Marko- 
vian and non-Markovian entanglement dynamics for dif- 
ferent choices of the parameters of the system and the 
bath (0J2/0J1, A/wj, 7/wi, hsT/cui, r). The following 
considerations are valid for separate baths as well as 
for a common environment. In the case of very small 
coupling (7 ~ 10 _3 wi) non-Markovian corrections are 
practically unobservable independently on the value of 
the system parameters up to relevant temperatures. To 
observe some appreciable deviations, we need to reach 
temperatures of the order of 10a>i. This result is under- 
stood considering that for small values of the system-bath 
coupling the role played by the bath itself is highly re- 
duced. Since non-Markovian corrections to the values of 
the coefficients are relevant only during the first stage of 
the evolution, it can be also predicted that if the initial 
state is robust enough (as, for instance, in the case of 
squeezing r = 2 discussed in the previous sections) , non- 
Mar kovianity will have marginal effects. It is in fact true 
that, in the case of 7 ~ 10~ 3 wi and T ~ IOcji, a correc- 
tion to the entanglement decay time can only be observed 
assuming low initial squeezing (r < 0.1). When 7 starts 
to increase, it is possible to deal with cases where some 
observable corrections due to non- Mar kovianity emerge 
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also in the low-temperature regime. In particular, the 
case of equal frequencies turns out to be the more sensi- 
tive. To give an example, in Fig. [7j we compared the 
two evolutions for two different values of r. Starting 
from r = 2, non-Markovian corrections are almost neg- 
ligible, giving corrections of few percents to the value of 
tF found within the Markovian approximation. In order 
to observe a noticeable difference, we must start with a 
factorized state (r = 0). As it can be observed in Fig. [7j 
in this case, while memoryless baths are able to support 
a very small quantity of entanglement, the initial kick 
given by non-Markovian coefficients enhances entangle- 
ment generation, its maximum amount is amplified of 
about one order of magnitude, and also its death time is 
increased. It is however worth noting that the absolute 
value reached is, in any case, very low. 





2 3 

W-|t 



FIG. 7: (Color online). Comparison between Markovian 
(orange) and non-Markovian (black) entanglement dynam- 
ics. The system parameters are ui = ui2, A = 0.2u}f, 
7 = 0.01 x 2/Wl, A = 20wi, and k B T = 10wi. The ini- 
tial squeezings are a) r = 2 and b) r — 0. 



V. DISCUSSION AND CONCLUSION 

During last years, a series of papers has investigated 
entanglement dynamics of coupled harmonic oscillators 
in contact with thermal environments. Nevertheless, a 
comprehensive analysis around the role played by the di- 
versity of the two oscillators was missing. 

We started our discussion with the case of separate 
baths, which seems to be more physically meaningful. 
As one should expect, after a transient, unless the tem- 
perature is very low, the initial two-mode squeezed state 
I^tms) becomes separable. Studying the decay time, 



we observed a scries of interesting features that we have 
qualitatively explained. First of all, by increasing u^/wi 
and by keeping unchanged all other parameters, entan- 
glement survival time increases as well. In the limit of 
UJ2/0J1 3> 1, fi+ and f2_ tend, respectively, to UJ2 and lo\. 
Since the asymptotic thermal state the system is reaching 
depends itself on the frequencies, the purity of the final 
state is enhanced together with the coherence time nec- 
essary to reach it. On the other hand, near u>2 = wi, en- 
tanglement survival time can be raised also by increasing 



|A|. When |A| 



fi_ is close to zero. Here, we argue 



that the behavior we observe is the result of two com- 
peting effects. If, on one side, the thermal equilibrium 
state is less pure because of the presence of a vanishing 
frequency, on the other side, the number of decohering 
channels goes to zero, since 0_ falls in a poorly popu- 
lated region of the spectral density, increasing in this way 
the time necessary to reach the thermal equilibrium. 

If the two oscillators share a common bath, it becomes 
possible the observation of asymptotic entanglement at 
relevant temperatures. One of the two eigenmodes of the 
system can result in fact decoupled from the environ- 
ment. As a consequence of this decoupling, the asymp- 
totic density matrix will result as the direct sum of the 
matrix of pure state and of a thermal state, correspond- 
ing to the mode interacting with the bath. When the 
two oscillators have identical frequencies and bath cou- 
plings, the mode X- gets decoupled independently of the 
value of A, giving rise to the plateau we observe in Figs. 
IH [SJ It is important to stress that this regime can be 
achieved not only when the oscillators are resonant, but 
also in presence of frequency detuning, through a proper 
engineering of the system-bath interaction. Conversely, 
if their bath couplings are different 71 7^ 72, frequencies 
uji and UJ2 can be found so that this special regime is 
achieved. 

We also studied the dynamical evolution of the degree 
of quantmnness of the system by means of the quan- 
tity d, introduced in Eq. (jT2"j) . The negativity of this 
quantity is a direct evidence of the negativity of the 
Glauber-Sudarshan quasi-probability distribution of the 
state. The agreement between the behavior of d and the 
entanglement is worth to be investigated, since a tight 
relation between them is not known. 

As for the Markovian versus non-Markovian debate, 
we compared the two behaviors in a wide range of pa- 
rameters. While in Ref. (27J the authors investigated 



what changes when the frequency of the oscillators goes 
outside from the spectral distribution of the bath in the 
high temperature regime, we tried to do an extensive 
analysis for LJi,ui2 *C A. Our results indicate that non- 
Markovian corrections can be observed, and are impor- 
tant, only for a reduced subset of initial conditions. To 
see them, we must put a relevant coupling between os- 
cillators and bath, u)\ should be close to uj 2 , and the 
initial state should be (almost) disentangled. Since the 
master equation has been obtained in the weak-coupling 
limit, the first of this assumptions seems at least ques- 
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tionable. If we release only one of them, Markovian and 
non-Markovian evolutions are almost identical. 
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and Xi(r) = a a xi + a i2 x 2 + Papi + (3 l2 p 2 , with 

an = cos 2 9 cos fl-T + sin 2 9 cos fl+r, (A4) 
sin 29 

a\2 = — - — (cosii + T — cosil-T), (A5J 

cos 2 6 . n sin 2 6 . n 
Pn = — — sinS2_rH — — — smS2+T, (Ao) 



fi 



+ 



sm29 sinfi+r sinf2_r 

2 ( n+ n~'' 



(A7) 



Appendix A: Derivation of master equation 

A system-bath model is usually described through the 
Hamiltonian H = Hs + Hb + V, where Hs refers to the 
system and Hb to the bath, and where V is the term 
containing the interaction. Generally, the coupling term 
can be written as V = J2k^k ® Bk, where Sk denote 
system operators and B^ bath operators. The master 
equation for the reduced density matrix to the second 
order in the coupling strength (for a general discussion 
see for instance j37j) is given by 



- \ dr^ {Cl, m (T)[SlS m (-T)p - S m {-T)pSi] 
Jo l,m 

+ C mil (-T)\pS m (-T)Sl-Sl P S m {-T)]}. (Al) 



and a 22 (9) = a u (7r/2 - 9), (3 22 (9) = Pu(ir/2 - 9), a 21 = 
ai2, and /3 2 i = fli2. The coefficients appearing in Eqs. 
(|A4HA7|) are defined as 



a 1 ( 2X 

9 = — arctan 



and 



± V 2 2 

The normal modes for the position are 

Q + = cos 9x 2 + sin 9x\ , 
Q- = cos #21 — sin 9x 2l 

while for the momentum 



(A8) 



(A9) 



(A10) 
(AH) 



Here, C; )TO (t) are the bath's correlation functions, de- 
fined by 



P + = cos 9p 2 + sin 9pi , 
P- = cos9pi — sm9p 2 . 



(A12) 
(A13) 



C hn (r) =Tr B [Bi(r)B m (0)Ro 



(A2) 



being Rq 



-0H b 



/Tre ^ Hb ) the equilibrium density 



matrix of the bath. In Eqs. (|AllA2p the tilde indicates 
the interaction picture taken with respect to Hq = H$ + 
H B - For instance, Bj(r) = e iH ° r B x eT iHaT . 

In the following, the correlation functions will ap- 
pear combined as Cf m (r) = Ci, m (r) + C mi ;(~r) and 
Cf m (r)=C l , m (r)-C m , l (-r). 



1. Separate baths 



In our model, whose Hamiltonian is H sep , we can iden- 



tify Si = £,A«X« and B2 = J2 k X^X^> and, cor- 
rcspondingly, Si = x\, S 2 = x 2 . Their expressions in the 
interaction picture are 



B, 



cos Q^'t + P k 



(») 



(A3) 



Notice that, for ui\ = uj 2i 9 = 7r/4. This implies an = 
a 22 and /?n = (3 22 . 

Since the two baths are identical and uncorrelated, we 
have Ci,i(t) = C 2 , 2 (r) = C(t) and C 1j2 (t) = C 2 ,i(t) = 
0. From the knowledge of the density of states J(fi) it is 
possible to obtain the explicit expression of C c (r) and 
C A (r): 



C°(r) = -i / dnj(Q) sin fir, 
Jo 



(A14) 



O 



poo 

C a {t) = / dfi J(S1) cos fir coth — — . (A15) 
Jo 2k B T 



The master equation reads 



dp 



[Hs,P\ - / dr ^2 {i x tA x t(~ T ),P}} 

J i=i 2 



+ [^[^(-rj.p]]^^}. 



C C (r) 
2 

(A16) 



By using expressions (| A4B A7|) . it can be written as in Eq. 
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(01), with 

4 = -* / dTa l3 C c {T), (A17) 
Jo 

Da = [ draijC^T), (A18) 
Jo 

Fij = I dTp l3 C A {r), (A19) 
Jo 

Tii = * / dTp t] C c {T). (A20) 
Jo 

The considerations about symmetry given in Sec. [[Jean 
be understood in terms of the explicit form of the mas- 
ter equation. The substitution A — > —A is equivalent to 
sending 6 to —6. Then, also ai2 and /3i 2 change their 
sign. But, since these coefficients multiply one operator 
(position or momentum) of the first oscillator and one 



operator of the second one, the canonical transformation 
U compensates the change of sign, and the evolution of 
p is unchanged. 

The explicit form of eij , Dij , Tij , Fij , in the Markovian 
case, can be calculated using the equalities 

lim / ri re -'("-"o)T = -uj ) + i P | (A21) 
where P denotes the Cauchy principal value, and 

+ oo 

™^-f-=2k B T J2 ~^T7a' ( A22 ) 
2kbJ ^ oj + vi. 

n— — oo ,b 

with v n = 2TmksT. 

From the master equation, a set of closed equations of 
motion for the average values of the second moments can 
be derived = 1,2): 



Jt 

dipiPj) 
dt 



d({xi,Pj}) 
dt 



1 

~2 



^(A + e\ 2 ) J2[({x k , Pi })0- - 8 ktj ) + ({x k , Pj })(l ~ Sk,i)] 
fe=i 

2 

){PiPj) - r i2 53[(Pfcj»*)(l - + (PkPj){l - S k ,i)] + Di 



fc=i 



2<p j p i ) - 2( Wj 2 + 4-)<^j> - 2 ( A + e ? 2 ) - <%) + {sis,)] 



(A23) 



(A24) 



(A25) 



2. Common bath 

In this case, C 1a (t) = C 2 , 2 (t) = Ci )2 (r) = C 24 (r). 
Then, the master equation can be written as has 



~dt 



[H s ,p}~ dT {[ X i>i x j{-T),p} 

170 «=1,2 



,C A (r) 



}■ 



2 

(A26) 



and, using Eqs. (|A17HA20|) . assumes the form given in 
Eq. © with Da = Du + D12 and similarly for e 2 ,T,F. 
The equations of motion for the second moments read 
(i,j = l,2) 



11 



d(xjXj) 
dt~ 

djpiPj) 
dt 



d ({xj,Pj}) 
dt 



(jPiXj + PjXi) 
2 

- ^2( u t + 4k)(( x kPj)S lk + (piX k )d jk ) 
fe=i 

1 2 1 



k=l 



(T u + Tjj)(j)iPj) -^rnKftp^l - S k j) + iP]Pk){l - hi)} 



k=l 



2(j3 lP2 ) -2(u? + e 2 jj )(x i x j ) 



2$> + %k){*i*k)Q- - - 2(A + 4 )(* 2 )(1 - 5 tJ ) 

fe=i 

2 

f 3 j ({sctjPi}) - X] ^({^i,Pfe})(l - *tfc) + Fa 



k=l 



(A27) 



(A28) 



(A29) 
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